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ABSTRACT 

We investigate how uncertainties in flux measurements affect the results from 
modified blackbody SED fits. We show that an inverse correlation between the 
dust temperature T and spectral index (3 naturally arises from least squares fits 
due to the uncertainties, even for sources with a single T and (3. Fitting SEDs 
to noisy fluxes solely in the Ray leigh- Jeans regime produces unreliable T and 13 
estimates. Thus, for long wavelength observations (A ^ 200 /xm), or for warm 
sources (T 60 K), it becomes difficult to distinguish sources with different 
temperatures. We assess the role of noise in recent observational results that 
indicate an inverse and continuously varying T — (3 relation. Though an inverse 
and continuous T — [3 correlation may be a physical property of dust in the ISM, 
we find that the observed inverse correlation may be primarily due to noise. 

Subject headings: dust - infrared:ISM - methods:miscellaneous 

1. Introduction 

Recent advances in infrared and sub-millimeter observations have enabled detailed in- 
vestigation of the properties of dust in a wide range of environments. Both ground and space 
based observatories, such as SCUBA, Bolocam, MAMBO, Spitzer, and ISO, as well as the 
balloon borne experiment PRONAOS, have revealed much about the nature of dust emis- 
sion. The near-future observatories Herschel and Planck will also be capable of detecting 
dust in a variety of environments. One conclusion drawn from spectral energy distribution 
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(SED) fits to fiux measurements by PRONAOS is that the dust emissivity spectral index 
j3 has an inverse correlati on with the dust t e mperature T, and t hat this correlation is well 
described as a hyperbola (jPupac et al.ll2003l ). iDesert et al. J2008h. using loriger wa velengths 
fluxes of galactic sources from the Archeops experiment, and I Yang &: Phillips! (120071 ). through 
observations of luminous infrared galaxies (LIRGS), also find a similar anti-correlation. 

In this short paper, we present an investigation of the effect of noise on the T and /3 
estimates from SED fits to fiux measurements. In recent work, we focused on starless cores, 
with T ^ 20 K, and showed that SED fits may result in inversely correlated and erroneous 
T and [3 estimates, due to line of sight temperature variations, noise, or both (Shetty et 
al., ApJ submitted, hereafter Paper I). Here, we consider warmer sources. We investigate 
how uncertainties in the fiuxes, coupled with the T — (3 degeneracy that aris es from the 



funct i onal form of the m odified blackbody spectrum describing dust emission flBlain et al. 



20031 : ISajina et al.ll2006l ). lead to a spurious inverse T — j3 correlation from SED fits. 



2. Modified Blackbody Assumption 



The common assumption abo ut the SED due to dust is that it is a blackbody modified 
by a power law in the frequency (IHildebrandl 1 19831 ) . For optically thin dust emission, the 
fiux density Sy takes the form 

S, = nB,{T)Ko(^^^\. (1) 

Here, Vl is the solid angle of the observing beam, B^{T) is the Planck function, and N is 
the dust column density. The term kq^u/uq)'^ is the frequency dependent opacity of the 
emit ting dust. Dust iii the interstel lar medium is usually characterized by 1 ^(3^2 



[e.g. iDraine fc Ledll984l : lMathislll990l . and references therein). Equation ([T]) only accurately 
characterizes sources with a single T (isothermal) and a single (3. 

With 4 or more fiuxes, a direct fit of equation ([1]) to the fiuxes can be employed to 
estimate three parameters: (3, T, and kqN (Paper I, and references therein). As discussed in 
Paper I, the fitting results are sensitive to noise and temperature variations along the line of 
sight. A single temperature fit to fiuxes near the peak of the SED may result in erroneous T 
and P estimates, due to line of sight temperature variations. For fiuxes in the Rayleigh- Jeans 
(R-J) tail, the best fit T may provide an accurate estimate for the "column temperature," or 
density weighted temperature. However, in the R-J tail, the fits are very sensitive to noise 
uncertainties. In this work, we consider isothermal sources, focusing on the effect of noise 
on T and (3 estimates from SED fits. 
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3. Effect of Noise on SED Fits 

In order to assess the effect of noise on tlie SED fits, we employ simple Monte Carlo 
experiments. We compute tlie fluxes at various wavelengths from sources with given tem- 
peratures and spectral indices, using equation ([1]). To incorporate the effect of noise, we 
then multiply each flux by a factor f+e, where e is a random value drawn from a Gaussian 
distribution with a chosen dispersion a (and mean of 0). We then perform a minimized Chi- 
squared fit of equation ([1]) to the noisy fluxes, and compare the resulting (3 and T estimates 
to those of the sources. 

Figure [T] shows the 75% and 50% probability contours of the best fit T and /? for three 
isothermal sources, with T=20, 60, and 100 K, and with (3=2. To obtain good statistics, 
we perform fits to numerous sets of data. The fluxes are sampled at 10 yum intervals, in 
the range A G 100-600 /im. The noise levels of the fluxes are 5% and 10%. For a 5% error 
in the fluxes, the mean (or median) of the T and (3 fits provide a good estimate of the 
source parameters for the 20 K source. Though an inverse T — (3 relationship is apparent, 
the spread about the source parameters is very small. However, for the warmer sources, 
the fits show a clear inverse T — (3 relationship; the mean or median T and (3 also do not 
correspond to the source values. We note that weighting the fluxes by signal-to-noise, or 
fitting in log(A)-log(S'i,) space, can improve the fits, but an inverse T — (3 relationship still 
emerges. In this work all the fluxes are weighted equally, in order for direct comparison with 
recent investigations in f|5l 

Also shown in Figure [T](b) are fits to fluxes from a 20 K source with /3=1.5. These 
fits illustrate that sources with different values of the spectral index, but with identical 
temperatures, would give vertically offset, but otherwise similar, points in a T — /5 plane. 

In practice, only a few fluxes are typically available from observations for fitting an 
SED. To simulate this more realistic scenario, we consider fits to sparsely sampled fluxes. 
Figure [2] shows the 75% and 50% likelihood contours of the best fit T and j3 from fits to five 
fluxes, with A=100, 200, 260, 360, and 580 /im. In comparison with Figure [H the spread in 
both (3 and T estimates is larger, especially for fits to fluxes from the 20 K source. The fits 
from the warmer 60 K and lOOK sources are less distinguishable from each other, compared 
with the analogous fits incorporating many more fluxes shown in Figure [H Further, fits to 
fluxes from the warmer sources are less likely to recover the source temperature: for the 60 
K source, and with cr=5%, only ~20% of the fits recover temperatures within 10% of the 
source temperature; for the 100 K source, only ~10% of the fits accurately recover the source 
temperature within 10%. 

Also shown in Figure [2] (as a dashed line) is the best fit hyperbola to PRONAOS and 
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IRAS data found by iDupac et al.l (120031 ) . The wavele ngths of the five flu xes considered in 
the fits in Figure [2] are the same as those considered by lDupac et ahl (120031 ) . After describing 
the degeneracy between (3 and T in the next section, we discuss imphcations of our analysis 
on recent observational interpretations in §3 



4. Degeneracy between the Spectral Index and the Temperature 



We have shown that modest levels of noise can result in grossly misleading T and (3 
estimates from a least squares fit of equation ([1])E The best fit is anti-correlated with 



Sajina et al. 


2006; 


Desert et al. 


2008) 



2003 



In the R-J part of the SED, fits to fiuxes from various warm sources may be difficult 
to distinguish (e.g. the 60 and 100 K fits in Figure [2]). The reason for the confusion in the 
resulting fits is illustrated in Figure [31 Figure M^a.) compares two SEDs, with T=60 K and 
29.3 K, and (3=2 and 3.4, respectively. The latter is a fit to five noisy fiuxes from the 60 
K source. The two SEDs are rather similar at large wavelengths, especially in the range 
A G 100-600 /im, as shown in Figure [3t^b). In this wavelength range, there are numerous 
combinations of T and (3 that are consistent with the observed SED, and therefore it is 
difficult to constrain them. Thus, fitting modified blackbody SEDs solely to fiuxes in the 
R-J regime produces unreliable T and f3 estimates. Moreover, the degeneracy between T 
and f3 is such that an overestimate in f3 results in an underestimate in T. 

For isothermal sources, fits to fiuxes near the peak of the SED provide more accurate 
estimates of the source parameters (Paper I). Thus, the 100-600 /im wavelength range con- 
sidered here can be shown to provide accurate estimates for source with T ^ 20 K. This 
is why the mean value of T and (3 recovers the source values for the cold 20 K source in 
Figures [1] and [21 but not so for the warmer sources. Shorter wavelength observations in the 
Wien regime may mitigate the effect of noise for warmer sources; but, at wavelengths A ^ 
100 /im, embedded sources as well a s transiently heate d very small grains (in cold regions) 
may contribute to the observed fiux (iLi fc Drainell200ll ). 



^Other techniques, such as Bayesian or hkelihood based methods, may provide better estimates of the 
source parameters. 
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5. Comparison with Recent Observations 



Pupae et al.l (120031 ) find an inverse correlation between tlie best fit T and (3 from 
PRONAOS data. In their analysis, they consider the degeneracy between T and (3 by 
performing statistical tests on noisy model fluxes. They find that noise is insufficient to 
account for the derived T — (3 trend, and suggest that an inverse T — (3 correlation is an 
intrinsic characteristic of dust in t he ISM. In our exte nsion of their work, we find a greater 
scatter in T and j3 (see Fig. 4 of iDupac et al.l (1200 ll )). That the best fits populate the T 
and j3 plane with a very similar shape to that produced simply by the presence of noise is 
strongly suggestive that noise is affecting these fitting results (see Figs [T] and [2]) . 



The T and (3 fits from lDupac et al.l (120031 ) show lower values of /3 at T ^ 20 K, compared 
with the distribution in Figures [1] and [21 As already discussed, comparing isothermal sources 
with different values of (3 would simply result in vertically offset points on the T — (3 plane. 
We thus consider whether sources with di fferent temperature s, but with /3 =1.5, are consisten t 
with PRONAOS data, focusing on Ml 7 JPupac et al.ll2002h and Orion JPupac ejallboOlh. 



which are the two sources that span the full range T G 20-80 K in Figure 3 of iDupac et al. 



(l2003f ). 



Figure H] shows numerous best fit T and j3 to five noisy fiuxes from 30 K, 40 K, and 50 
K isothermal sources, wi th (3=1. b. For the warmer sources, t he fits are similar t o many of 
the points in Figure 4 of lPupac etaP J2002h at T > 20 K. As lPupac et~aD J2002h describe, 
at some positions the IRAS 100 /xm fiuxes were not included in the fit. Figure H] also shows 
such fits. Though this simple model cannot account for the points at low T in Figure 4 of 
Pupae et al.l (|2002| ). many of the points at T ^ 20 K are very similar. Thus, a model with 
constant (3 ~ 1.5, for dust at temperatures ^ 30 K, may account for many of the PRONAOS 
derived T and (3 points (at T ^ 20 K). Such temperatures may be realistic for much of the 
dust in M17 (e.g. IColdsmith et al.lEOT . 



To verify whether any substantial differences can be easily identified in cases where j3 
is constant or variable, we perform the Monte-Carlo experiments described in §3] assuming 
models with different (3 — T dependencies. The f irst model ha s /?=1 .5 for all T; the second 
has /3=l/(0.4+0.008T), which is the fit found by lPupac et al.l (120031 ). We construct a series 
of sources with T G 15-80 K, in increments of 5 K, and fit equation ([1]) to numerous sets of 
noisy fiuxes at 100, 200, 260, 360, and 580 /im. Figures El^a)-(b) show the results for fits 
to fiuxes wit h nois e levels of 5%, for both models, along with the best fit hyperbola from 
Pupae et all J2003h . 



Figures [5](a)-(b) show that the fit parameters from the lowest temperature sources at 
15, 20, and 25 K are distinct from the warmer sources. As already discussed, at these 
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temperatures, given the wavelengths of the fluxes considered, the peaks of the emergent 
SEDs are well sampled, so the fits are not very sensitive to noise. 

The main difference between the constant-/^ and the inverse- hyperbolic /3(T) models 
is that the there is less scatter in the fit (3 from the model where (3 depends on T. This 
difference occurs because the sources with low temperatures have higher spectral indices 
(and vice versa); thus, this inverse-hyperbolic (3{T) function reduces the scatter about [3 for 
a given T relative to the scatter from the constant (3 case. 

The extent of the scatter in Figure[5]is dependent on the chosen level of noise, the form of 
/3(T), as well as the temperatures of the sources considered in this test (and the wavelengths of 
the sampled fluxes). For example, for the range T G 40 — 50 K in the constant-/? model (green 



points in Fig. [S](a)), the fits are well described by the best fit hyperbola from iDupac et al. 



(120031 ). Further, simply increasing the noise level in the model where [3 varies inversely with 
T (Fig[5](b)) could produce a.T — (3 scatter plot that is similar to the constant 13 model with 
the same temperature range (Fig. [5]^a)). The primary effect of an intrinsic inverse T — (3 
correlation is to reduce the spread in the T — j3 plane due to noise in the flux measurements. 

Despite the differences in the scatter between the models in Figure [5]^a)-(b), there is 
a striking similarity in the shape of the T — 13 distribution: even though the models have 
different forms in l3{T), fits to the fluxes result in a similar inverse T — (3 correlation. Due 
to this similarity, it is very difficult to distinguish an intrinsic T — (3 anti-correlation from 
the artificial anti-correlation arising due to noise. Alternative estimates of the actual range 
in source temperatures, such as molecular line observations or radiative transfer modeling, 
would be necessary to accurately account for the spread in T and (3 due to the intrinsic 
degeneracy between the parameters in least squares SED fitting. 

One method to quantify the contribution of measuremen t uncertainties to de rived cor- 



relations is the use of correlation coefficients (e.g. iKellyl 120071 ). iDupac et al.l (1200 ll ) (in §4.2) 
compare the correlation coefficient of their best fits from Orion observations to a model with 
no correlation between T and j3. The model has numerous uncorrelated T and (3 pairs, from 
which noisy fluxes are constructed; fits to those noisy fluxes result in a correlation coefficient 
between T and 13 of -0.4. The discrepancy between this value and the value from the ob- 
servations, -0.92, i ndicates that a uniform and uncorrelated distribution in T and (3 can be 



ruled out (see also lYang &: Phillips! 120071 ). Such a test, however, has not ruled out scenarios 



where (3 is constant, the temperature range is limited, or both. 

The correlation coefficients from T — (3 fits to noisy fluxes from any isothermal model are 
^ -0.90, which are remarkably similar to the value from the fits to Orion observations. In a 
case where there is a small range in T, such as the green points in Figure [5](a) (T G 40-50 
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K), we compute a correlation coefficient of -0.8, also similar to the value found from the 
observations. These statistical tests suggest that a constant /3, with a limited range in the 
source temperature, cannot be ruled out in interpreting the PRONAOS Orion observations. 

If the P-ROA^-A 06* obs ervations detected sources with T ^ 20 K, and given the relative 



flux uncertainties of ~5% (IDupac et al.ll2001uPajot et al.ll2006l ). then the lack of points with 
T ^ 20 and /9 ^ 1.2 is suggestive that /? for such cold sources cannot be a constant value of 
1.5. One possibility is that P=2 for the coldest sources, but P = 1.5 for sources with T ^ 30 
K. Similar to an inverse hyperbolic function in such a step function would result in a 

T — P distribution with less scatter compared to a scenario where /3 is constant. Figure [5]^c) 
shows the fits from this step function in /3(T). This model illustrates that other functional 
forms of /3(T) may produce a T — /3 distribution similar to that found by the PRONAOS 
observations. 

The observationally inferred T—f3 correlation is certainly consistent with, e.g. an inverse 
intrinsic hyperbolic P(T) function. Yet, our analysis suggests that the intrinsic functional 
form of /3(T) cannot be reliably derived from a fit to the T and /3 values obtained from a least- 
squares SED fit to the observed fluxes. Other functional forms of /3(T), in combination with 
the noise levels, may also be consistent with the observations. The agreement between the 
T — (3 correlation coefficients from fits to PRONAOS Ohon observations and our isothermal, 
single P models suggest that T and /5 may not be correlated. Other statistical tests are 
needed to quantify the T — /5 degeneracy. 

A more complete understanding of the distribution in the source T and P is necessary 
to accurately infer any physi cal T — co r relatio n. The scatter in the best fit parameters 



about the hyperbola in the iDupac et al.l (120031 ) study is likely a consequence of various 
factors. Since the wavelengths of the observed fluxes are in the R-J part of the emergent 
SED from w arm sources, noise d oes significantly bias the fitting. Additionally, the dataset 
presented in iDupac et al.l (120031 ) contains numerous sources, with different spectral indices 
and temperatures, between sources and perhaps even within individual sources; these issues 
may all contribute to the scatter in the distribution of T and p. 



6. Summary 

We have shown that a natural consequence of fitting modified blackbody SEDs to flux 
measurements, with modest noise uncertainties, is an inverse correlation between the tem- 
perature T and spectral index p. Such an inverse correlation arises even for isothermal 
sources with a constant p. Least squares fits to fluxes from the Rayleigh- Jeans regime of the 
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emergent SED are very sensitive to noise uncertainties. Consequently, various isothermal 
sources with different Ts but identical /3s may be indistinguishable through simple SED fits. 

We find that the spurious T — (3 ant i- correlation due to noise, from numerous isothermal 
sources with a limited range in T and /3, is similar to the observationally inferred anti- 
correlation (§5). A hyperbolic T — (3 relation may indeed be a physical characteristic of dust, 
but the observational results show similarities with other explanations, such as a single (3 
for isothermal sources with T ^ 30 K, or a step function where (3^2 for T ^ 20 K and 
f3 «i 1.5 for T ^ 30 K. Thus, least squares SED fits to measured fluxes may not be able to 
reveal any intrinsic correlation between f5 and T. 
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Fig. 1. — The 75% and 50% probability contours for the best fit T and /? to noisy fluxes, in 
the wavelength range 100-600 /xm (sampled at 10 /xm increments), from 20 K, 60 K, and 100 
K isothermal sources with 13 — 2 (marked by crosses). Gaussian distributed noise is included 
in each flux, with (a) a—b% and (b) a—lQi%. In (b), best flt T and (3 from a 20 K source 
with is also shown. 
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Fig. 2. — The 75% and 50% probability contours for the best fit T and /3 to noisy fiuxes, 
with (a) o"=5% and (b) cr=10%. The wavelengths of the sampled fiuxes are A=100, 200, 260, 
360, and 580 /zm, from three isotherrn al sources (as in Fig. [1]). The dashed line is the best 
fit hyperbola from lDupac et al.l (120031 ). 
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Fig. 3. — Actual and fit SEDs from a 60 K isothermal source. The boundary of the shaded 
region is the dust SED from a 60 K source. The solid curve shows a fit to the few noisy 
fluxes (crosses) at A=100, 200, 260, 360, and 580 //m. The errors in those fluxes are 2.1%, 
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Fig. 4. — Best fit (3 and T to noisy fluxes (with a=5%) from isotliermal sources witli T g30- 
50 K. Colored points sliow flts to fluxes witfi A=100, 200, 260, 360, and 580 /im. Black points 
sliow flts to fluxes from a 30 and a 40 K source excluding the 100 /xm flux. The horizontal 
line indicates the spectral index of th e source, /3=1.5. The dashed line shows the best flt to 
data presented by lDupac et al.l (120031 ). 
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Fig. 5. — Best fit (3 and T to noisy fluxes (with cr=5%) from isotfiermal sources witfi a range 
of temperatures, for three sets of models, along with the best fit from iDupac et al.l (120031 ) 
(dashed line), (a) f3=1.5; (b) /5=l/(0.4+0.008T); and (c) a step function in (3, with (3=1.5 
for T > 20 K and p=2.0 for T=20 K. The fiuxes have A=100, 200, 260, 360, and 580 ^im. 



